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ABSTRACT 


A  Massey  diffuse  convolutional  decoder  was  simulated  in  the  IBM 
7030  computer  raid  used  to  decode  binary  digital  error  patterns  meas¬ 
ured  at  2400  bits/sec  on  an  HF  radio  circuit  and  on  a  multiple  link 
data  circuit  dominated  by  a  tropospheric  scatter  path.  It  vas  found 
that  the  simulation  performed  at  least  as  veil  as  previous  cyclic 
code  simulations  on  HF  radio  and  better  than  cyclic  codes  on  the 
tropospheric  path. 
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SECTION  I 
INTRODUCTION 


The  wide-scale  application  of  high-speed  digital  computers  for 
information  processing  in  data  management  and  process  control  systems 
has  resulted  in  an  increased  demand  for  digital  communications.  This 
demand  has  been  met  by  the  development  of  digital  data  transmission 
devices  for  the  existing  analog  communications  plant. 

Due  to  the  inherent  redundancy  of  most  analog  communication, 
the  quality  of  the  transmission  media  is  not  particularly  critical 
to  the  attainment  of  satisfactory  information  exchange.  However,  in 
the  case  of  digital  communication  an  unacceptable  occurrence  of 
errors  may  be  introduced  in  the  transmission  process  if  suitable 
redundancy  is  not  incorporated  into  the  transmitted  data  stream. 

The  attainment  of  acceptable  digital  error  performance  using  cable, 
microwave  and  communication  satellite  links  has  been  achieved  prima¬ 
rily  by  careful  circuit  engineering  with  provision  of  adequate  signal 
to  noise  margins.  On  media  such  as  high-frequency  radio  and  tropo- 
scatter  the  same  quality  of  link  engineering  cannot  always  be  achieved. 
Consequently,  the  incorporation  of  error  control  has  been  of  particular 
interest  with  emphasis  on  applications  to  digital  data  transmission 
on  poor  quality  links.  This  error  control  is  effected  through  the 
addition  of  controlled  redundancy  to  the  transmitted  data  stream  in 
the  form  of  the  binary  digital  coding  technique  which  permits  the 
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detection  and/or  correction  of  errors  by  the  error  control  device 
installed  at  the  receive  terminal. 

Among  the  coding  techniques  that  have  been  developed  for  error 
control  purposes  is  the  Massey  diffuse  convolutional  code  [l].  In 
this  paper,  the  results  of  a  computer  simulation  of  the  performance 
capabilities  of  this  coding  technique  against  measured  real  channel 
error  patterns  is  investigated.  The  use  of  a  valid  computer  simula¬ 
tion  permits  a  controlled  (repeatable)  investigation  of  the  perform¬ 
ance  of  error  control  in  a  manner  which  would  not  be  possible  in 
field  environment  tests. 

The  Massey  Diffuse  Code 

Hie  Massey  diffuse  code  [l]  evaluated  in  this  paper  is  a  l/2 
rate  code  (i.e.  the  code  supplies  one  redundant  bit  per  source  in¬ 
formation  bit. )  The  encoder  accepts  from  the  source  information 
bits  (i)  at  a  clock  interval  (t).  Since  the  error  data  is  for  2400 
bit/sec  channels  and  the  code  rate  is  l/2,  the  clock  interval  (t) 
corresponds  to  1200  bits/sec.  Redundant  (parity)  bits  are  calculated 
according  to  the  relationship: 

^t  -  +  *t-x  +  i~t-2x.  +  it-(  3x+l)  (1) 

where  i  represents  an  information  bit 
p  represents  a  parity  bit 

x  is  some  specific  number  of  clock  intervals  and  addi¬ 
tion  is  mod  2  with  no  carry. 
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The  pair  of  bits  and  i^  are  then  transmitted  over  the  channel. 

The  received  bits  (both  information  and  parity)  are  entered 
into  the  decoder  (Fig.  l).  In  the  decoder  syndromes  ( S)  are  formed 
in  accordance  vith  the  following  relations: 


Sx+1  '  ell  +  eix+l  *  epxtl 
S2x+1  +  S3x+1  “  eil  *  ei3x+l  +  ep2x+l  +  ep3x+2 

S3x+2  =  611  +  eix+2  +  ei2x+2  +  ei3x+2  +  ep3x+2 

.  f  1  if  the  bit  was  in  error 

where  ei,  or  ep^.  =  ■<-.. 

t  t  L  0  otherwise 


(2) 


The  fact  that  the  error  bits  are  used  is  seen  by  examining  the 
decoder  structure  and  equation  (l).  If  equation  (l)  is  satisfied, 
the  output  of  the  initial  information  and  parity  addition  will  be 
zero.  The  decoding  rule  is  to  change  i^  if  the  sum  of  equations 
(2)  is  greater  than  2.  This  criterion  will  be  correct  if  not  more 
than  two  of  the  eleven  bits  in  equations  (2)  are  in  error. 

The  encoder  when  supplied  with  an  information  message  of  all 
binary  O's  will  calculate  the  parity  to  be  all  binary  0Ts.  This 
allows  for  the  simplification  of  the  simulation  in  that  if  all  zero 
information  is  encoded  and  the  bit-by-bit  channel  errors  are  added 
to  the  encoder  output  mod  2  with  no  carry,  the  result  is  the  error 
patterns.  Thus,  it  is  only  necessary  to  provide  the  error  patterns 
as  input  to  a  simulated  decoder  and  record  the  errors  which  are  output 
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Figure  1  Massey  Diffuse  Convolutional  Decoder 


from  the  decoder. 


The  decoder  while  correcting  errors  also  introduces  delay  into 
the  communication  system  to  which  it  is  applied.  This  delay  is  a 
function  of  the  parameter  x  ( identified  as  B  in  the  simulation 
program)  and  is  presented  in  Table  I. 

Description  of  the  Data 

The  error  data,  have  been  described  in  detail  previously  [2,  3> 

4  ]  and  will  be  briefly  reviewed  here.  In  October  of  1965,  field 
tests  were  conducted  on  NKD  HF  circuits  between  Antigua  and  Ascension 
Islands  to  collect  HF  digital  error  patterns  [2,  3]*  Tests  were  con¬ 
ducted  for  10  minute  periods  (runs)  at  2400  bits/sec  using  a  Kineplex 
TE-21 6  phase- shift-keyed  l6  tone  modem.  The  procedure  was  to  trans¬ 
mit  a  52  bit  test  message  from  Antigua  to  Ascension  and  back  to 
Antigua  on  a  looped  basis.  At  Antigua  the  received  message  was 
compared  with  a  suitably  delayed  replica  of  the  52-bit  test  message, 
and  bits  which  were  not  in  agreement  were  declared  in  error.  The 
results  were  recorded  on  magnetic  tape  and  converted  to  IBM  seven- 
track  format  magnetic  tape  for  use  in  the  IBM  7030  computer.  The 
other  data  was  collected  [4]  in  the  spring  of  1966  on  a  link  (Fig.  2) 
characterized  as  having  three  types  of  transmission  media:  wireline, 
microwave  and  troposcatter.  The  troposcatter  is  multiple-hop,  and 
the  wireline  consists  of  the  interconnection  of  numerous  leased 
telephone  wirelines.  The  dominant  sub-path  is  a  troposcatter  hop  of 
583  miles.  During  the  test  Sebit  24b  vestigial  sideband  AM  modems 
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Figure  2  Full  Duplex  Communication  Circuit  Configuration 


Table  I 


DELAY  DUE  TO  DECODING 


B=X 

Delay 

6 

.015 

l4o 

•  350 

144 

.360 

284 

•  710 

288 

.720 

572 

1.^3 

576 

1.44 

1196 

2.99 

1200 

3-oo 
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were  vised  at  2^00  bits/ sec  in  a  h  KHz  channel  which  was  5M  multiplexed 
with  other  channels  for  transmission  on  the  dominant  sub-path.  The 
test  procedure  was  the  same  as  for  the  HF  tests  except  that  the  trans¬ 
mission  was  on  a  one  way  basis  and  runs  were  90  minutes  long. 

Analysis  of  the  HF  data  [2,  3]  indicates  that  errors  occur  in 
bursts  and  are  in  no  way  random.  A  typical  burst  is  from  two  to 
four  thousand  bits  long  and  has  an  error  density  of  5%  Bursts  are 
separated  by  intervals  which  have  a  low  error  density  and  vary  in 
length  from  one-tenth  the  length  of  the  burst  they  follow  to  one- 
million  times  the  burst  length.  Additionally,  periodic  errors  of 
spacing  32  bits  occur  due  to  tone  interference  on  individual  tones 
of  the  TE-216  modem.  The  data  collected  on  the  multiple-hop 
path  consisted  of  bursts  that  varied  in  length  from  a  few  hundred 
bits  to  tens  of  thousands  of  bits  where  for  the  most  part  the  density 
of  errors  was  75$*  Hie  intervals  between  bursts  were  error- free  and 
varied  from  one-tenth  the  length  of  the  burst  they  followed  to  one 
million  times  the  length. 
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SECTION  II 
SIMULATION  RESULTS 


Each  of  the  HF  channel  test  runs  listed  in  Table  II  vas  processed 
through  the  computer  simulation  of  the  Massey  diffuse  decoder.  The 
output  error  patterns  were  recorded  on  magnetic  tape  for  further 
statistical  processing.  The  gross  results  of  performance  on  the  HF 
data  will  be  presented  in  terms  of  output  average  error  rate  versus 
B  in  Table  II  and  as  improvement  factor  in  Figures  3  through  7. 


vhere  Improvement  Factor  (F)  =  Channel  Average  ErrorRate  £e) 

Output  Info.  Average  Error  Rate  ( 


V 


Values  of  and  Pc  are  presented  in  Fortran  floating  point  nota¬ 

tion.  The  decoder  gives  from  one  to  three  orders  of  magnitude 
improvement  with  100$  error  correction  on  numerous  runs.  Further, 
the  improvement  increases  with  B  (i.e.  as  delay  is  increased  per¬ 
formance  increases).  Hie  values  of  B  were  selected  to  give  a 
range  of  delay  and  were  adjusted  such  that  bits  on  the  same  tone 
were  not  encoded  together.  This  eliminates  the  possibility  that 
the  periodic  errors  will  affect  an  information  bit  and  all  of  its 
associated  parity  bits  in  the  channel  and  thereby  cause  decoding 
errors.  In  no  case  was  the  output  average  error  rate  (P  )  greater 
than  the  channel  average  error  rate  (Pg) .  Hie  decoder  did  not 
degrade  the  communication  channel. 

The  results  for  the  tropospheric  data  are  presented  in  Table  III 
and  Figures  8  through  12.  Since  this  data  was  collected  on  a  channel 


(3) 
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Table  II 


OUTPUT  ERROR  RATE  AS  A  FUNCTION  OF  B  IN  THE  HF  CHANNEL 

Input 

Error  Output  Error  Rate 


Run  # 

Rate 

B=6 

B=l40 

333 

3.  5E-5 

0 

0 

324 

7.6E-5 

0 

0 

291 

7.8E-5 

0 

0 

318 

1.2E-4 

2-3E-5 

0 

321 

1. 3E-4 

0 

0 

345 

6.0E-4 

l. 5E-4 

1.4E-5 

354 

6.IE-4 

7.7E-5 

8-9E-6 

248 

6.2E-4 

1-7E-4 

0 

258 

7.6E-4 

l.OE-4 

6.  5E-5 

315 

1.7E-3 

3-45E-4 

2.4E-4 

342 

2.8E-3 

l.OE-4 

4.6E-5 

264 

2.9E-3 

3.0E-4 

1.  4e-4 

309 

4. 3E-3 

1.8E-4 

9-5E-5 

261 

9-2E-3 

2. IE-3 

1.  6e-3 

312 

9-3E-3 

2.2E-3 

9.6e-4 

249 

l.OE-2 

1.2E-3 

9-5E-4 

339 

1.2E-2 

2.7E-3 

1.4e-3 

B=284  B=572  B*1196 

OOO 
0  0  0 

0  0  0 

0  0  0 

0  0  0 

1-7E-5  6. 5E-6  l.OE-5 

1.  IE- 5  8.9E-6  6.0E-6 

0  0  0 

7-7E-5  7-3E-5  5- 

2.2E-5  9-6E-5  5- 

4.6e-5  4.2E-5  7- 

1.7E-4  4.0E-5  2. 

7-2E-5  6.8E-5  2. 

l.  5E-3  l.  5E-3  8. 

l.OE-3  9.6E-4  6. 

9.0E-4  l.OE-3  9- 

l.  3E-3  1-2E-3  l. 
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Figure  3  Simulation  Test  Results,  B  =  6 
Improvement  Factor  (F)  Versus  Channel  Average  Bit  Error-Rate  (Pp) 


11 


If) 

<D 

rO 

OJ 

I 

< 


HF  CHANNEL 


Figure  4  Simulation  Test  Results,  B  =  l40 
Improvement  Factor  (F)  Versus  Channel  Average  Bit  Error- Rate  (P  ) 
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Figure  5  Simulation  Test  Results,  B  =  284 
Improvement  Factor  (F)  Versus  Channel  Average  Bit  Error-Rate  (P  ) 
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Figure  6  Simulation  Test  Results ,  B  =*  572 
Improvement  Factor  (F)  Versus  Channel  Average  Bit  Error-Rate  ( Pe ) 
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Figure  7  Simulation  Test  Results,  B  =  1196 
Improvement  Factor  (F)  Versus  Channel  Average  Bit  Error-Rate  (Pg) 
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Table  I IT 


OUTPUT  ERROR  RATE  AS  A  FUNCTION  OF  B  IN  THE 
TROPOSPHERIC  CHANNEL 


Run  # 

Input 

Error 

Rate 

B=6 

6o 

1.2E-5 

9-7E-6 

64 

2.0E-5 

3.6E-6 

6i 

2.4E-5 

1. 3E-5 

56 

2-7E-5 

1. 3E-5 

49 

3- IE- 5 

l.  8e- 5 

42 

4.  4e-5 

3-9E-5 

4i 

5-  3E-5 

4.  6e-5 

33 

6.2E-5 

6.  IE- 5 

43 

6.8E-5 

6.7E-5 

57 

7.2E-5 

2.6E-5 

30 

7-9E-5 

6.4e-5 

51 

8.OE-5 

4.3E-5 

70 

8.OE-5 

6.8E-5 

69 

8. 4e- 5 

2.5E-5 

46 

9.8E-5 

5-8E-5 

62 

i.oe-4 

9-3E-5 

63 

l.OE-4 

l.OE-4 

53 

l.IE-4 

6.OE-5 

38 

i.ie-4 

6.9E-5 

4o 

1.4e-4 

l.IE-4 

52 

1. 5E-4 

3.  IE- 5 

71 

1. 6e-4 

l.IE-4 

55 

1-7E-4 

7-9E-5 

59 

1-7E-4 

1.7E-4 

39 

2.0E-4 

l.IE-4 

37 

2.0E-4 

l.OE-4 

66 

2.1E-4 

2.  3E-4 

32 

2.2E-4 

2. 5E-4 

31 

2.6E-4 

2. 3E-4 

68 

2.8E-4 

l.IE-4 

36 

3-OE-4 

1.  3E-4 

65 

3-7E-4 

3-7E-4 

48 

4.0E-4 

4. 3E-4 

44 

7-OE-4 

5-9E-4 

34 

8. 3E-4 

8.7E-4 

67 

9-OE-4 

l.OE-3 

35 

9-3E-4 

l.OE-3 

Output  Error  Rate 


B=l44 

B=288 

B=576 

i.4e-6 

0 

0 

0 

0 

0 

4.  5E-6 

0 

0 

1.  5E-4 

0 

0 

9.IE-6 

5.9E-6 

2.6e-6 

0 

0 

0 

l.IE-6 

2.2E-6 

1.6E-6 

6.2E-5 

4. 5E-5 

2.8E-5 

3-5E-5 

2.9E-5 

1.2E-5 

6.  6e-6 

4.3E-6 

2.  IE-6 

2.8E-5 

l-  3E-5 

2.OE-5 

l.  6e- 5 

1-3E-5 

1.3E-5 

2.OE-5 

i.oe-6 

4.0E-5 

7-9E-6 

8.8e-7 

0 

2.5E-5 

2.6E-5 

1.4E-5 

7- IE- 5 

3-6E-5 

8.3E-6 

6. oe- 5 

2.8E-5 

1.2E-5 

5-4E-5 

5-2E-5 

3.2E-5 

1. 6e- 5 

1.  IE- 5 

6.2E-7 

5.2E-5 

3.4E-5 

2.9E-5 

1.  6e- 5 

l.  6e- 5 

l.IE-5 

5.8E-5 

3-8E-5 

2. 4e- 5 

6.OE-5 

4.2E-5 

2.9E-5 

l.  5E-4 

1. 6E-4 

l. 6e-4 

3.4E-5 

3*4e- 5 

1. 3E-5 

5-7E-5 

5-9E-5 

2.6E-5 

2.  3E-4 

l.IE-4 

4.8E-5 

2.7E-4 

2.8E-4 

2.6E-4 

2.  5E-4 

2.0E-4 

i.ie-4 

3-3E-5 

3.0E-5 

2.3E-5 

9-6E-5 

l.OE-4 

7  •  4e- 5 

1.9E-4 

1.4E-4 

9 . 4e- 5 

5.  6e-4 

4.0E-4 

9.IE-4 

3. 3E-4 

2-9E-4 

2.2E-4 

8. 3E-4 

7-9E-4 

7.0E-4 

7-2E-4 

6.2E-4 

5.7E-4 

1.1E-3 

l.OE-3 

8.9E-4 
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B=l200 

0 

0 

0 

0 

0 

0 

0 

2.2E-6 

9-7E-7 

o 

1-3E-5 

o 

2.0E-5 

0 

2.7E-6 

1.2E-5 

5.2E-7 

9-2E-6 

0 

5-3E-7 

2.9E-7 

1- 9E-5 
1.1E-5 
8.4e-5 
4.8e-6 

o 

3-7E-5 

9-5E-5 

8.0E-5 

2- 3E-5 
6.9E-5 
7-5E-5 
5- IE-4 
2.  IE- 4 

3-  6e-4 
5-  IE- 5 
6.2E-4 


Table  III  (Continued) 


OUTPUT  ERROR  RATE  AS  A  FUNCTION  OF  B  IN  THE 
TROPOSPHERIC  CHANNEL 


Input 

Error  Output  Error  Rate 


Run# 

Rate 

b=6 

B=l44 

B=288 

B=576 

B=1200 

47 

1.1E-3 

1.2E-3 

1.2E-3 

l.IE-3 

l.OE-3 

8.9E-4 

50 

1.2E-3 

l. 3E-3 

1.2E-3 

1.2E-3 

1. 5E-3 

1.4E-3 

58 

2.2E-3 

2.4e-3 

2.4E-3 

2.2E-3 

1.9E-3 

1.4E-3 

45 

2.4E-3 

2.5E-3 

4.7E-3 

4.7E-3 

4.7E-3 

4.7E-3 

54 

4.0E-3 

4.1E-3 

4.1E-3 

4.2E-3 

3.8E-3 

3.0E-3 
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Figure  8  Simulation  Test  Results,  B  ^  6 
Improvement  Factor  (F)  Versus  Channel  Average  Bit  Error-Rate  (Pg) 
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Figure  9  Simulation  Test  Results,  B  =  1^^ 

Improvement  Factor  (f)  Versus  Channel  Average  Bit  Error-Rate  (Pg) 
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Figure  10  Simulation  Test  Results,  B  =  288 
Improvement  Factor  (F)  Versus  Channel  Average  Bit  Error-Rate  (Pe) 
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Figure  11  Simulation  Test  Results,  B  =  57° 
Improvement  Factor  (F)  Versus  Channel  Average  Bit  Error-Rate  (P  ) 
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Figure  12  Simulation  Test  Results ,  B  =  1200 
Improvement  Factor  (F)  Versus  Channel  Average  Bit  Error-Rate  (P^) 
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with  no  periodic  errors,  there  was  no  need  to  specially  adjust  the 
values  of  B.  It  is  obvious  from  the  scatter  diagrams  that  the  value 
b=6  was  insufficient  for  the  long  dense  bursts  in  that  improvement 
was  less  than  one  order  of  magnitude  in  error  rate  and  on  many  runs 
P^  >  P^  (alternately  F  <  l).  However,  as  B  increases  performance 
likewise  increases,  and  when  B  reaches  1200  (3  sec  delay)  only 
three  runs  show  F  <  1  while  the  other  runs  have  up  to  three  orders 
of  magnitude  improvement  with  twelve  runs  completely  corrected.  It 
is  obvious  from  the  scatter  diagram  that  this  trend  would  continue 
if  B  were  further  increased. 

Cumulative  Performance  Curves 

A  common  method  for  presenting  the  performance  of  a  communica¬ 
tion  system  in  which  errors  occur  is  to  present  the  percent  of  channel 
data  intervals  for  which  the  error  rate  is  better  than  a  particular 
value  of  error  rate.  This  can  be  done  for  a  coded  system,  and  the  curves 
can  then  be  compared  to  the  uncoded  channel.  In  order  to  properly 
compare  decoded  to  uncoded  performance,  the  decoded  interval  length 
must  be  that  length  which  contains  the  information  portion  of  the 
corresponding  channel  data  interval. 

Ulus, 

Decoded  Channel  Interval  Length  _ 

Uncoded  Channel  Interval  Length  ~  o  a  e  (4) 

These  results  are  presented  in  Figures  13  and  l4  for  the  HF  and 
Tropospheric  data  respectively. 
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Figure  13 


Cumulative  Performance  of  Massey  Code  in  HF  Channel 


2k 


Figure  14  Cumulative  Performance  of  Massey  Code  in  Tropospheric  Channel 


SECTION  III 


CONVOLUTIONAL  CODE  ERROR  PROPAGATION 
It  is  known  from  the  work  of  Massey  that  convolutional  codes,  when 
they  fail,  propagate  errors.  Further ,  at  some  stage  in  the  operation 
of  the  decoder  it  is  no  longer  necessary  for  channel  errors  to  occur  in 
order  to  have  the  propagation  continue.  This  is  called  infinite  propa¬ 
gation.  An  advantage  of  a  computer  simulation  is  that  it  is  possible 
to  tag  the  channel  errors  and  follow  them  through  the  decoder.  By  do¬ 
ing  this,  it  is  possible  to  determine  which  errors  were  not  corrected, 
what  caused  any  additional  output  errors,  and  whether  or  not  infinite 
propagation  occurred.  In  Tables  IV  and  V  the  percent  of  output  errors 
that  are  uncorrected  input  errors  for  the  HF  data  and  tropospheric 
data  respectively  is  presented.  All  of  the  remaining  errors  are  gener¬ 
ated  errors  with  no  direct  counterpart  in  the  channel  error  patterns. 

The  generated  errors  were  compared  to  the  source  error  table,  and  it 
was  found  that  in  all  cases  the  generated  error  was  within  a  distance 
B  of  the  source  channel  error.  Thus,  error  generation  occurred  when 
a  burst  affected  a  coded  unit  of  information  bits,  and  there 

was  no  spillover  beyond  the  code  length.  An  alternative  way  of  stat¬ 
ing  this  is  that  if  errors  occurred  in  three  of  the  four  coded  informa¬ 
tion  bits ,  an  error  was  generated  into  the  fourth  bit  which  due  to  the 
diffusion  of  the  code  occurred  B  bits  later  in  time.  There  was  no 
infinite  propagation  observed.  In  fact  a  proof  [l]  exists  that  infinite 
propagation  cannot  occur  for  this  particular  code  structure. 
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Table  IV 


°!o  OF  OUTPUT  ERRORS  THAT  ARE  UNCORRECTED  INPUT  ERRORS 

HF  DATA  RUNS 


Run 


Number 

b=6 

B=l4o 

B=284 

B=572 

B=1I 

333 

— 

_ 

_ 

_ 

_  - 

324 

-- 

-- 

-- 

-- 

-- 

291 

-- 

— 

-- 

318 

50 

-- 

-- 

-- 

— 

321 

— 

— 

— 

— 

3^5 

77 

44 

75 

67 

60 

354 

54 

33 

75 

33 

100 

248 

60 

-- 

-- 

-- 

-- 

258 

53 

52 

54 

54 

42 

315 

52 

57 

75 

85 

57 

342 

54 

60 

63 

46 

57 

264 

53 

64 

59 

64 

54 

309 

64 

59 

63 

67 

4o 

26l 

55 

56 

55 

58 

56 

312 

53 

56 

58 

58 

64 

249 

58 

53 

55 

56 

53 

339 

48 

52 

51 

52 

52 
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Table  V 


1°  OF  OUTPUT  ERRORS  THAT  ARE  UNCORRECTED  INPUT  ERRORS 
TROPO  DATA  RUNS 


Run 


Number 

b=6 

B-144 

60 

64 

75 

64 

69 

6l 

66 

63 

56 

66 

50 

49 

59 

62 

42 

67 

— 

4i 

79 

50 

33 

72 

68 

43 

77 

67 

57 

65 

51 

30 

65 

67 

51 

60 

65 

70 

71 

62 

69 

69 

58 

46 

64 

58 

62 

59 

57 

63 

65 

60 

53 

62 

62 

38 

73 

51 

4o 

67 

58 

52 

53 

54 

71 

57 

54 

55 

60 

55 

59 

78 

78 

39 

72 

61 

37 

67 

6l 

66 

76 

6 ^ 

32 

74 

70 

31 

69 

63 

68 

71 

60 

36 

64 

58 

65 

66 

63 

48 

74 

11 

44 

60 

58 

34 

67 

59 

67 

75 

65 

35 

73 

63 

47 

62 

60 

50 

74 

77 

58 

67 

63 

45 

84 

82 

54 

67 

64 

B=288  B=576  B=1200 


71 

57 

75 

100 

_  _ 

71 

78 

67 

65 

64 

100 

78 

54 

-- 

52 

55 

59 

72 

65 

— 

80 

51 

57 

75 

-  - 

66 

62 

78 

57 

67 

55 

61 

59 

67 

65 

62 

53 

57 

50 

— 

54 

53 

62 

63 

100 

55 

55 

55 

65 

67 

37 

73 

68 

59 

65 

63 

69 

60 

5 8 

61 

64 

66 

70 

66 

66 

60 

52 

60 

52  . 

60 

DO 

54 

59 

58 

62 

63 

60 

74 

68 

59 

60 

58 

57 

59 

56 

60 

65 

53 

57 

66 

66 

60 

64 

63 

56 

77 

70 

64 

64 

63 

64 

79 

74 

65 

61 

62 

60 
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SECTION  IV 


C  CMP  ARISON  OF  THE  CONVOLUTIONAL  TECHNIQUE  TO  CYCLIC  CODIM} 

Since  there  is  a  great  deal  of  knowledge  on  the  performance  of 
the  (24,  12,  3)  cyclic  Golay  code  on  the  measured  channels  [2,  5]  an 
idea  of  the  significance  of  the  results  for  the  Massey  code  can 
be  obtained  by  finding  the  performance  of  the  Golay  code  on  the  test 
data  for  the  same  amount  of  delay  used  with  the  Massey  code.  The 
(24,  12,  3)  code  can  correct  3  random  errors  in  24  bits  at  code  rate 
l/2.  The  code  will  be  interleaved  by  coding  together  bits  spaced  in 
bits  apart  to  obtain  a  code  word.  The  total  delay  introduced  by  this 
technique  is  24  m/l200  seconds.  If  m  =  150  then  the  delay  is 
3  seconds  and  a  total  interleaved  code  block  contains  24  x  150  =  3^00 
bits.  Generated  errors  can  occur  if  there  are  more  than  three  errors 
in  a  code  word  and  will  occur  within  that  code  word  with  no  overflow 
from  one  code  block  to  another.  Thus,  whatever  generated  errors  do 
occur  will  be  constructed  in  the  same  way  that  the  generated  errors 
of  the  Massey  code,  in  fact,  occurred.  The  performance  of  the  (24, 
12,  3)  code  was  found  by  simulation  [6]  for  a  delay  of  three  seconds 
(m  =  150).  The  performance  of  the  code  on  the  HF  data  was  found 
previously  [5]  and  will  not  be  reproduced  here.  Comparison  to  the 
previously  presented  results  indicates  that  the  Massey  diffuse  con¬ 
volutional  code  performes  as  well  as  the  Golay  code.  The  results 
found  with  the  tropospheric  data  and  presented  in  Table  VI  and 
Figure  15  are  far  different.  For  this  channel  the  performance  of 
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Table  VI 


PERFORMANCE  OF  GOLAY  CODE  ON 
TROPOSPHERIC  CHANNEL 


Run  # 

Output 

BER 

F 

Run  § 

60 

2.14E-5 

0.5 

71 

64 

3.00E-7 

66.6 

55 

ol 

3-00E-7 

80 

59 

56 

1.40E-5 

1.9 

39 

49 

3-OOE-7 

103 

37 

42 

3.00E-7 

146 

66 

4i 

2.15E-6 

24.6 

32 

33 

1.94E-4 

•31 

31 

43 

3-OOE-4 

22  6 

68 

57 

3.00E-7 

24o 

36 

30 

1.40E-5 

5-6 

65 

51 

4.00E-5 

2 

48 

70 

1.15E-6 

69.5 

44 

69 

5-OOE-6 

16.8 

34 

46 

2.15E-6 

45.5 

67 

62 

4.00E-6 

25 

35 

63 

1.20E-5 

8.3 

47 

53 

6.20E-7 

177 

50 

38 

3-  30E-5 

3.3 

58 

4o 

1.10E-5 

12.  Y 

^5 

52 

1. 15E-6 

130 

54 

Output 

BER 

F 

2.00E-5 

8 

3-OOE-7 

5  66 

2-97E-4 

0-5 

6. 51E-4 

•  3 

1.08E-3 

.18 

1.25E-4 

1.6 

1.75E-4 

1.2 

8.50E-5 

3-0 

3-00E-7 

993 

7-OOE-5 

4.2 

1. 50E-4 

2.4 

5-86e-4 

.68 

1.61E-4 

4. 3 

2.40E-5 

34.5 

7-43E-4 

1.2 

5-35E-4 

1-7 

I.27E-3 

.86 

2-39E-3 

0-5 

3.OOE-6 

73-3 

4.95E-3 

.48 

1-39E-3 

2.8 
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TROPOSPHERIC  CHANNEL 


Figure  15  Simulation  Test  Results  of  Interleaved  Cyclic  Code  with 

Three  Seconds  Delay 

Improvement  Factor  (F)  Versus  Channel  Average  Bit  Error-Rate  (P  ) 
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the  Massey  diffuse  convolutional  code  is  substantially  superior  to 
that  of  the  Golay.  The  Massey  code  achieved  a  zero  output  error 
rate  on  12  of  the  42  runs  while  the  Golay  code  did  not  achieve  a 
zero  output  error  rate  on  any  of  the  runs.  The  Massey  code  demon¬ 
strated  overall  error  generation  (F  <  l)  on  3  of  the  42  runs  while 
the  Golay  code  did  this  on  9  of  the  42  runs.  In  Figure  1 6  a  compari¬ 
son  is  made  between  the  cumulative  performance  of  the  Massey  diffuse 
code  with  B  =  1200  (3  sec.  delay)  and  the  interleaved  cyclic  Golay 
code  with  the  same  delay. 
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Fifure  16  Comparative  Performance  of  Massey  and 
Cyclic  Codes  for  Equal  Delay 
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SECTION  V 


CONCLUSIONS 

'"ho  objective  of  this  paper  has  been  to  obtain  the  performance 
of  a  Massey  diffuse  convolutional  code  on  actual  measured  HF  and 
tropospheric  error  patterns*  Based  on  the  data,  three  general  con¬ 
clusions  can  be  reached. 

1)  Improvement  increases  with  delay  with  no  upper  bound  on 
performance  within  the  range  of  delays  considered 
(Delay  £  3  sec. ) . 

2)  Error  generation  occurs  over  a  limited  distance  from 
the  excessive  burst  of  errors  causing  the  generation. 

3)  a)  On  the  HF  data  the  Massey  code  performance  is  as  good 

as  that  of  the  Golay  code  , 

b)  On  the  Tropospheric  data  the  Massey  code  performance 
is  superior  to  that  of  the  Golay  code- 
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